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Abstract 



Recently it has been pointed out that diagonal Pade approximants to 
truncated perturbative series in gauge theories have the remarkable property 
of being independent of the choice of the renormalization scale as long as the 
gauge coupling parameter a(p 2 ) is taken to evolve according to the one-loop 
renormalization group equation - i.e., in the large- j3$ approximation. In this 
letter we propose and describe an improvement to the method of diagonal 
Pade approximants. The improved method results in approximants which are 
independent of the chosen renormalization scale when a(p 2 ) evolves at any 
chosen (in principle arbitrary) loop-level. 
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Pade approximants (PA's), either diagonal or nondiagonal, when applied to any trun- 
cated perturbative series (TPS), possess by construction the same formal accuracy^ as the 
original TPS. This means that, when expanding the applied PA in a power series of the 
expansion parameter of the TPS, we reproduce all the coefficients of the TPS. PA's, being 
rational fractions, additionally act as a kind of analytical continuation of the TPS and thus 
often represent a substantial improvement of the results deduced directly from the TPS. 
This had been the main motivation for applying PA's to TPS's in gauge theories such as 

QCD i.g 

Recently it has been noted [Q that the diagonal Pade approximants (dPA's), when 
applied to TPS's in gauge theories, have in addition the remarkable property of being 
independent of the choice of the renormalization scale (RScl) if the gauge coupling pa- 
rameter a(p 2 ) is taken to evolve according to its one-loop renormalization group equation: 
a{p 2 ) = a(q 2 )/[l + /3 \n(p 2 /q 2 )a(q 2 )}. This is the direct consequence of the mathematical 
property of dPA's that they are invariant under the homographic transformations of the 
argument (Ref. [lj, Part I): z^azj{\ + bz). Since a full observable (formally an infinite 
perturbation series) must be RScl-independent, the mentioned property of dPA's suggests 
that the dPA for an available TPS of an observable in a gauge theory (approximately) sums 
up a substantial set of diagrams and thus represents a very reasonable resummation method 
there. The authors of || further investigated this dPA method and showed that the re- 
summed diagrams represent systematic approximations to the Neubert's || concept of the 
distribution of momentum flow through a bubble-dressed gluon propagator. The authors 
of [H pointed out the need for further improvements of the method, in particular to obtain 
RScl-independence beyond the large-/?o (one- loop flow) limit. 

We present here an algorithm which improves the dPA method in this sense - by con- 
structing approximants which are RScl-independent at any chosen loop-level of evolution of 
a(p 2 ) and which, when expanded in power series of a, give the same formal accuracy as the 
TPS's to which they are applied. 

A generic observable S in a gauge theory (e.g., QED or QCD) can in general be redefined 
so that it has the following form as a formal perturbation series: 

S = a{q 2 )f{q 2 ) = a{q 2 ) [l + ri (q 2 )a(q 2 ) + r 2 {q 2 )a 2 {q 2 ) + ■■■ + r n {q 2 )a n {q 2 ) + ■ • 

Here, a(q 2 )=a(q 2 )/7r and q 2 is a chosen renormalization scale (RScl). The full series ([!]) is 
of course independent of q 2 . In practice we have only a limited number of coefficients rj(q 2 ) 
available (i = l,. . . ,n), i.e., we know only a truncated perturbation series (TPS) 

S n {q 2 ) = a{q 2 )f n \q 2 ) = a{q 2 ) [l + ri (q 2 )a(q 2 ) + r 2 {q 2 )a 2 {q 2 ) + ■■■ + r n (q 2 )a n (q 2 )} . (2) 

This TPS explicitly depends on the RScl q 2 - changing q 2 changes the value of S n in general 
by a term ~ a n+2 . The RScl-independence of ([D and the RScl-dependence of (|2|) suggest 



1 We refer to Ref. (Part I) for the basic theory of PA's. 

2 For a review of methods of dealing with power expansions in quantum field theory, see Ref. 
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that an approximant that has the same formal accuracy and is RScl-invariant is a good 
candidate to be closer than (g) to the full answer ([[]). 

We will now construct such approximants. The gauge coupling parameter a(p 2 ) = a(p 2 ) / it 
evolves according to the perturbative renormalization group equation (RGE) 



da(p' 
din p 



, -£/V' + V), 

j=0 



(3) 



where f3j are constants if particle threshold effects are ignored. They are positive in QCD 
and negative in QED. Only a limited number of these perturbative coefficients (5j (/3 ,. . . ,^3) 
are known in QCD (cf. 0, in MS scheme) and in QED (cf. ||, in MS, MOM and in on-shell 
schemes). Hence, in practice the RGE (|3|) will always be truncated at some level. We now 
introduce the ratio of gauge coupling parameters at two different renormalization scales p 2 
and q 2 



k(a q , u) 



a(p 2 
a(q 2 



where: a q = a(q ) , u = \n(p /q). 



(4) 



Formally expanding this function in powers of u = \n(p jq ) results in the following series: 



k(a q ,u) = 1 + £ u 3 kj(a q ) , where: kj( 



1 &> 



We note that kj(a q )^a^(q 2 ) since RGE (^|) gives the connection 

k,(a q ) = (-lypla^q 2 ) + O (a 3+1 (q 2 )) , k (a q ) = 1 



(5) 



(6) 



and the terms of higher orders can also be explicitly obtained from RGE ([|). At this point 
we rearrange the formal power series ([!]) for S/a(q 2 ) into a related series in kj(a(q 2 )) 



S = a{q 2 )f{q 2 ) = a(q 2 ) 



3=1 



(7) 



,rj(q 2 ) 



We note that the coefficient fj{q 2 ) depends solely on the first j coefficients ri(q 2 ). 
of the original series (|l|), as implied by relations (||). In addition, fj(q 2 ) depends on the first 
j coefficients fio,- ■ ■ of RGE (f|). We then define the corresponding formal series ^(q 2 ), 
which is in powers of (— /3o(i(q 2 )) 



a{q 2 )T{q 2 ) = a(q 2 



3=1 



(8) 



We construct for aiq^Fiq 2 ) the diagonal Pade approximants (dPA's) with argument a(q 2 ) 



a(q 2 )[M-l/M]Aq 2 ) = a(q 2 



A/-1 



1 + E a m (q> m (q 2 " 



m=l 



M 



1 + Y: b n (q 2 )a n (q 2 ) 



n=l 



-1 



a(q 2 )F(q 2 ) = a(q 2 )[M - l/M] r {q 2 ) + O (a 2M+1 (q 2 )) 



(9) 
(10) 
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The above dPA depends only on the first (2M— 1) coefficients fj(q 2 ) (j = 1,. . . ,2M — 1), 
due to the standard requirement (|10|) . Since, as mentioned earlier, the coefficient fj(q 2 ) is 
a unique function of only ri(q 2 ),. . . ,rj(q 2 ), the above dPA depends only on the first 2M — 1 
coefficients ri(q 2 ),. . . ,t 2 m-i( ( ? 2 ) of the original series (P, i.e., it is uniquely determined once 
we have the TPS (0) with n = 2M — 1 available. Unless we have an exceptional situation 
when the denominator in the dPA (||) has multiple zeros as polynomial of a(q 2 ), we can 
uniquely decompose this dPA into a sum of simple fractions 

Here, [— Ui{q 2 )j3o\~ l are the M zeros of the denominator of the dPA @ which is regarded 
as a polynomial of a{q 2 )- The above expression (|11]) is a weighed sum of the one-loop- 
evolved gauge coupling parameters^ a(pf), with generally complex scales p 2 determined by 
the relation Ui(q 2 ) = ln(p 2 /q 2 ), and with weights on. The approximant that we are looking for 
is then obtained by replacing in the above weighed sum the one-loop-evolved gauge coupling 
parameters with those which evolve according to the full RGE (|3[) 

M 

a(q 2 )Gf- 1/A \q 2 )^a(q 2 )Y,a t k(a(q 2 ),u i ) . (12) 

i=i 

The functions k(a(q 2 ),Ui) appearing here are defined via (f|) as ratios of gauge coupling 
parameters a at the RScl q 2 and the new scales p 2 = q 2 exp[ui(q 2 )} 

k (a(q 2 ),u^ = a(p 2 )/a(q 2 ) where: \n(p 2 /q 2 ) = Ui(q 2 ) . (13) 

We stress that a(pf) = a(p 2 )/ii — . . ,M) are the gauge coupling parameters evolved from 
the RScl q 2 to p 2 by the RGE (|J) whose loop-level precision (i.e., the number of included 
coefficients (3j) can be chosen as high as possible^, independently of the number n = 2M— 1 of 
the coefficients of the available TPS S n (q 2 ) of Eq. (0). By (|13"D, the obtained approximants 
(|j"2l) can be written in a somewhat more transparent form 

M 

a(q 2 )G^~ 1 ^ M \q 2 ) = ^ a^p 2 ) where: p^ = g 2 exp Ui(q 2 ) . (14) 

i=i 



Function k(a,u) in (|l2l), which depends on two (in general complex) arguments, can be 
called the kernel of the above approximant. We will call the above approximant the modified 
diagonal Baker-Gammel approximant (modified dBGA) with kernel k, since there exists a 
certain (but limited) similarity with the diagonal Baker-Gammel approximants as defined 
in Ref. |l| (Part II, Sec. 1.2). We emphasize again that this modified dBGA of order 



3 Evolved from the RScl q 2 to pf, by the one-loop ("large-/?o") version of RGE (|3|). 

4 In QED and QCD, this would mean inclusion of /3q,. . . ,/?3 since these coefficients are now available 
in certain schemes - see earlier discussion. 
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2M—1, for observable S of Eq. ([J), is uniquely determined by the first 2M — 1 coefficients 
ri(g 2 ),. . . ,r 2 M-i(<? 2 ) of the perturbative series ([!]), i.e., by the TPS SW-i of @ (n = 2M— 1), 
and by the independently chosen loop-level for the evolution of a(p 2 ) according to RGE 

Having the value a(q 2 ), we should stress that exact values a(p 2 ) can be obtained only if 
we evolve the RGE (|3]) numerically from u = ln(g 2 /g 2 ) = to u = ln(p 2 /q 2 ) = Ui(q 2 ). This 
numerical integration should be performed with additional care when £tj(g 2 )'s are complex. 
Approximate values for a(pf) can be obtained by perturbatively expanding the solution 
a(p 2 )/a(q 2 ) in powers of a(q 2 ), but such expansion would be reasonable only if the chosen 
RScl q 2 is not far away from the scales p 2 . In any case, the resulting approximant ( [HQ does 
not represent an analytical formula once we go beyond the large-/?o approximation. 

It can be shown that this modified dBGA ( |14|) of order 2M—1, for observable S of Eq. ([[]), 
fulfils the two requirements that we wanted to achieve: 

1. It has the same formal accuracy as the TPS .SW-i^ 2 ) of (0): 

S = a{q 2 )Gf- l/M \q 2 ) + O (a 2M+ \q 2 )) . (15) 



2. It is fully invariant under the change of the renormalization scale q 2 . In fact, the 
weights (ki and the scales p 2 = q 2 exp[ttj(g 2 )] are separately independent of the chosen 
renormalization scale q 2 . 

The formal proofs of these two statements will be given in a longer paper |J. Furthermore, 
also the discussion of similarities and differences between the presented modified dBGA's 
(Fj]) and the usual dBGA's of Ref. will be given in that longer paper. 

Within the presented algorithm, the case of one-loop evolution of a(p 2 ) (the large-/3o 
approximation) means: k(a, u) = l/(l+/3 u a), and kj(a) —/3q(— a)- 7 , when using notation of 
Eqs. (§)-(§)■ Therefore, in the one-loop case, expansion (|7]) for f (q 2 ) = S / a(q 2 ) and (H) for 
jF(g 2 ) are identical. The modified dBGA flT2"| ) is in this case reduced to the usual dPA ([□]). 

One may worry what happens when the parameters on and Ui{q 2 ) in the modified dBGA 
(Q)-(irj) are not simultaneously real. In that case, the modified dBGA could be complex. 
Since relation (|15"D and the RScl-invariance are valid for the entire modified dBGA's, they 
are valid separately for their real and imaginary parts. The observable S is real, so we then 
just take the real part of expression (|TJ). Since a(q 2 ) and 5* are real, relation ( |I5D implies 
that the imaginary part of the modified dBGA a(q 2 )G^~ 1 ^ M \q 2 ) must be ~a 2M+1 (g 2 ) or 
even less. 

It may be helpful not to remain at this rather abstract level, but to write more explicit 
formulas for dBGA's (|T2|)- (|T4]) in the practically interesting cases of M = l and M = 2. 

In the case M = 1 (n = 2M — 1 = 1) the method gives the same result as the effective 
charge (ECH) method and the Brodsky-Lepage-Mackenzie (BLM) method [12 in the 
large-/5 approximation 



Si(q 2 ) = a(q 2 )fV(q 2 ) = a(q 2 ) 1 + r 1 (g 2 )a(g 2 ) 



a(q 2 )G l f /1] {q 2 ) = a{Q 2 ) where: Q 2 



q 2 exp 



-ri(q 2 )/Po 



(16) 



Here, a(p 2 ) = a(p 2 )/ir evolves according to RGE @ where the chosen loop-level is arbitrary. 
It is straightforward to check directly that a(Q 2 ) is RScl-invariant and that (|l^) is satisfied. 



5 



In the case M = 2 (n = 2M— 1 = 3) , parameters of the dBG A (|T|)- (|14]) can also be obtained 
in a straightforward, although algebraically more involved, manner. RGE @ implies 

fci(a) = — /3 a — Pia 2 — /3 2 a 3 — • • • , (17) 
fo(a) = + (5/2)/3 0/ 9ia 3 + . . . , k 3 (a) = -^a 3 - . . . . (18) 



Here we use the short-hand notation a = a(q 2 ) = a(q 2 )/ix. Inverting relations (p~7|) ([T8|) 
yields expressions for a,a 2 and a 3 in terms of fci(a),/c 2 (a) and k^(a). We insert then these 
expressions into the truncated series S3 (q 2 ) of Eq. (|J) (if it is available) and thus obtain the 
first three coefficients of the rearranged truncated series for S^q 2 ) of the form (|?|) 



fs(q 2 ) 



/ 2 (g 2 ) = -§r 1 (g 2 ) + lr 2 (g 2 ), (19) 
Po Po Po 

(-ii + l) r '<« 2 ) + i^« 2 )-^ 2 '- < 2 °' 



With these coefficients, we form the truncated series for a(g 2 )J-"(g 2 ) (|8|), and the dPA 
a(g 2 )[l/2]jr(g 2 ) (|TTJ) in the form ([□]). We then obtain expressions for parameters Ui(q 2 ) 
and &i (i — l,M) for the case of M = 2 



u 2 ,i 



where: det 



(is - /1/2) ± v/det 



[2(/ 2 -A 2 )] ' (U2-U1 



\U2 - h) 1 „ , 01 x 

«i = 7- , OLi — 1 — cm . (21) 



/ 3 + /i(2/ 2 -3/ 2 ) ■ + 4(/ 2 - A 2 ) 3 . (22) 



The plus sign in corresponds to m 2 ((? 2 ). For simplicity, we omitted notation of the RScl- 
dependence in the coefficients fi(q 2 ) and v>i(q 2 ). Of course, expressions (|19|)-(f20|) should 
be inserted into (pT|) - (p2|) in order to obtain these parameters explicitly in terms of the 
original coefficients rj(g 2 ) (i = 1,2,3). When we insert these obtained parameters into the 
dBGA expression ([TJ]) (M = 2), we get the RScl-invariant approximation to S3 {q 2 ), with 
the RScl-invariant parameters <Sj and p 2 = q 2 exp[£tj(g 2 )] (« = 1,2) explicitly dependent on the 
original coefficients r^q 2 ) {i = 1,2,3) and on the RGE coefficients f3j (j = 0,1,2) of Eq. ([|). 
Although the parameters obtained above for the case of M = 2 contain dependence on the 
first three RGE beta-coefficients (flofli,^), we should emphasize that the evolution of the 
gauge coupling parameters a(p 2 ) appearingf^in the dBGA (|TJ]) can be governed by the RGE 
(0) with a higher chosen loop-level accuracy, e.g., by inclusion of (3% there. On the other 
hand, at least the first three coefficients (/3 ,/3i,/3 2 ) should be taken into account in the 
RGE evolution of a(p 2 ) from the RScl q 2 to pf since these three coefficients appear in the 
parameters <Sj and p 2 = q 2 exp[ttj(g 2 )] (i = 1,2). 

Concerning relations (|2~T|)-(|2"2"|) (for M—2), we can distinguish several cases 



When (/ 2 — ff) >0, then: Ui,ai are real (i = 1,2), u\^ui and 0<aj<l. 



RGE evolution from RScl q 2 to the (possibly complex) scale pf is meant here. 
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When (/ 2 -/ 2 )<0 and \f 3 + Wff -3f 2 )\>2^ (ft- f 2 )*, then: u h a t are real (i = l,2) 
and u\^u 2 . 



• When (/ 2 -A 2 ) <0 and | / 3 + A (2 A 2 - 3/ 2 ) | < 2^{f 2 -f 2 )\ then: u t are complex, a, 
generally complex = 1,2) , and U\^u 2 . 

• When (/ 2 -A 2 ) <0 and \f 3 + f x (2f 2 -3f 2 )\ = ^(f!-f 2 f [or when (/ a - / 2 ) = and 

fi]i then: the system of equations for «j and on is not solvable, i.e., form fllT] ) is 
not valid, the dPA (|^) has a multiple (double) pole. 

• When (f 2 — fi)—0 and f3 = ff, then: are real (i = 1,2) and u\=u 2 = f\. 

Even when the parameters <Sj and are complex [i.e., when 'det' in ( pT|) is negative], it can 



be checked directly from ( |2l| ) that 0:2 = (<5i)* and £t2 = (wi)*, and thus that the approximant 
is again real [note: a{{p\)*)=a{p\)*}. 

In QCD, presently available results of perturbative calculations contain, for various ob- 
servables 5* and in specific renormalization schemes, the coefficients ri(q 2 ) and r 2 (q 2 ) of (HI), 
but not yet r^{q 2 ). Hence, the described algorithm still cannot be applied for M = 2 for 
QCD observables, due to the fact that r 2 M-i(<? 2 ) = r 3 (q 2 ) are not yet known. This, however, 
is in stark contrast with some QED observables for which perturbative coefficients r 3 (g 2 ) 
have already been obtained. 

The presented algorithm, although being a clear improvement of the method of dPA's 
for perturbative series in gauge theories, still has several deficiencies. One of them is that 
the obtained approximants probably cannot discern in QCD, on the basis of a given TPS, 
the nonperturbative behavior originating from (ultraviolet und infrared) renormalons - see 
arguments in Refs. 0, ||]- |J for the case of the usual (d)PA's with which the presented 
modified dBGA approximants are closely related. Another deficiency is that the algorithm 
can be applied only in the cases when the available TPS S n (q 2 ) of (0) has odd n = 2M— 1, i.e., 
it cannot be applied in the case n = 2 (which is at present the case of many QCD observables). 
This is so because the algorithm heavily relies on the decomposition (|11|) which is valid only 
for diagonalPA's. This problematic restriction is present also in the large-/3o limit, i.e., in the 
usual dPA approach. Another problematic point is that the dBGA (P^)-(fT^) becomes in the 
cases of M>2 also explicitly renormalization scheme (RSch) dependent, because parameters 
Ui and oii appearing in such dBGA's also involve some of the f3j (j > 2) RGE coefficients 
which are, in contrast to (3q and (3\, RSch-dependent.Q This contrasts with some other 
approaches. For example, in the approach of the principle of minimal sensitivity (PMS) ||T0| . 
both the RScl- and RSch-independence of the approximant are achieved via a local method, 
while in the present approach the RScl-invariance is ensured via a more global method and 
RSch-invariance (i.e., independence of /?2,/?3,. . .) is not ensured at all. It would definitely 
be instructive to compare the efficiency of the presented method (for M = 2, for the time 
being in QED only) with the PMS method, as well as with other methods, among them: the 



effective charge (ECH) method |T|, the Brodsky-Lepage-Mackenzie (BLM) approach and 



6 When all effects of /3j's (j >2) are neglected, changing the renormalization scheme is equivalent 
to changing the renormalization scale. 
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its extensions ]T2[ — [14j] and 0, a new approach [TjJ based partly on the ECH methods, and 



another new approach |T6| using a method of analytic continuation. It may be also useful to 
investigate which classes of Feynman diagrams the presented algorithm approximately sums 
up,[] in analogy with the work || for the case of the dPA approach. 



Abbreviations used frequently in the article: 

(d)BGA - (diagonal) Baker-Gammel approximant; (d)PA - (diagonal) Pade approximant; 
RSch - renormalization scheme; RScl - renormalization scale; TPS - truncated perturbation 
series. 
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Under the term "approximately" we understand here a systematic approximation which would 
converge to the sum of mentioned diagrams when M^oo. 
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